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PROGRAM  SUMMARY 

This  report  summarizes  work  done  in  an  ongoing  program  under  contract 
RAD  80007  for  the  periods  October  1980  through  September  1981  (Project 
Order  RADC  81-10)  and  October  1981  through  September  1982  (Project  Order 
RADC  82-09).  The  objective  is  to  gain  information  on  the  interaction  of 
ionizing  radiations  with  various  communications  materials.  Thus,  the 
intent  is  that  work  performed  under  this  contract  will  increase  our 
understanding  of  the  basic  mechanisms  involved  when  photons,  electrons, 
or  other  ionizing  radiations  interact  with  matter  in  its  condensed 
phases.  Emphasis  is  being  placed  on  studying  the  interactions  occurring 
in  those  materials  which  are  currently  used,  or  which  are  being  considered 
for  use,  in  electronics  systems.  Knowledge  of  quantities  such  as  the 
cross  sections  for  the  various  types  of  electron  interactions  in  these 
materials  enables  predictions  to  be  made  concerning  the  performance  of 
electronics  systems  in  various  situations,  such  as  exposure  to  the 
results  of  a  nuclear  explosion.  Our  approach  to  this  general  problem  is 
both  theoretical  and  experimental.  In  our  experimental  program,  we 
continue  to  make  observations  which  can  be  related  to  photon  and  electron 
interactions  with  solids.  The  theoretical  portions  of  the  program 
include  predictive  calculations  for  a  variety  of  interaction  phenomena 
as  well  as  interpretative  studies  in  support  of  our  experimental  efforts. 
Results  from  the  basic  theoretical  and  experimental  programs  are  employed 
in  applied  calculations  with  direct  relevance  to  the  performance  of 
electronics  systems  including  microsystems  produced  by  VLSI  technology. 

In  studies  carried  out  during  this  2-year  period,  we  have  developed 
various  models  for  calculating  the  energy  loss  per  unit  length,  or 
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stopping  power,  and  the  inelastic  mean  free  path  for  electrons  in  a 
material.  We  have  included  both  organic  materials  such  as  polyethylene 
and  inorganic  materials  such  as  Si02-  In  each  case,  the  calculated 
values  of  stopping  power  and  mean  free  path  have  been  compared  with 
experimental  values  where  available  and,  in  the  case  of  mean  free  paths 
for  organic  materials,  with  those  predicted  from  our  universal  formula. 

In  addition  the  various  models  for  stopping  power  and  mean  free  path 
have  been  compared  with  each  other.  Measurements  have  been  made  of 
electron  attenuation  lengths  (=  electron  mean  free  paths)  as  a  function 
of  incident  electron  energy  for  electrons  in  amorphous  carbon  films,  as 
a  direct  check  on  the  values  predicted  by  our  universal  formula  for 
electron  mean  free  paths  in  solid  organic  insulators.  A  theory  has  been 
developed  to  describe  charged-particle  track  structure  in  nonmetal  lie 
solids,  as  the  track  evolves  in  time  and  space.  Work  has  continued  in 
the  field  of  microdosimetry  with  the  calculation  of  the  energy  densities 
deposited  by  high  energy  photons  as  a  function  of  depth  in  traversing  a 
Si 02  slab  sandwiched  between  two  slabs  of  Si. 

The  studies  reported  here  were  conducted  by  E.  T.  Arakawa,  J.  C.  Ashley, 
R.  N.  Hamm,  R.  H.  Ritchie,  J.  E.  Turner,  M.  W.  Williams,  and  H.  A.  Wright 
of  the  Health  and  Safety  Research  Division,  ORNL,  V.  E.  Anderson,  Computer 
Sciences  Division,  ORNL,  and  T.  A.  Callcott,  C.  Martin,  and  L.  R.  Painter, 
Department  of  Physics  and  Astronomy,  University  of  Tennessee,  Knoxville, 
Tennessee. 

Finally,  we  gratefully  acknowledge  the  continuing  advice  and  support 
of  0.  N.  Bradford  of  the  Rome  Air  Development  Center. 
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I.  ENERGY  LOSSES  AND  MEAN  FREE  PATHS  OF  ELECTRONS 

Our  studies  of  the  interaction  of  low-energy  electrons  with  solids 
began  over  a  decade  ago  and  since  then  results  from  this  work  have  been 
used  in  a  wide  variety  of  areas.  It  is  particularly  significant  to  note 
that  our  calculations  of  stopping  powers  of  solids  for  electrons  with 
energies  slO  keV  form  an  important  part  of  a  new  report  prepared  for 
submission  to  the  International  Commission  on  Radiation  Units  and 
Measurements  (see  M.  J.  Berger  and  S.  M.  Seltzer,  "Stopping  Powers  and 
Ranges  of  Electrons  and  Positrons,"  NBSIR  82-2550,  August  1982).  During 
this  reporting  period,  these  studies  have  continued.  The  results  are 
described  below. 

A.  Silicon  Dioxide 

Earlier  studies  of  electron  interactions  with  Si02  (RADC-TR-76-125, 
April  1976)  have  been  re-evaluated  and  refined  for  use  in  the  microdosi- 
metric  studies  undertaken  in  this  reporting  period.  New  theoretical 
results  for  stopping  powers  and  mean  free  paths,  and  comparisons  with 
experimental  data,  are  described  in  detail  in  the  following  references: 

1.  "Interaction  of  Low-Energy  Electrons  with  Silicon  Dioxide," 

J.  C.  Ashley  and  V.  E.  Anderson,  J.  Electron  Spectrosc.  Relat. 

Phenom.  24,  127-148  (1981). 

2.  "Energy  Losses  and  Mean  Free  Paths  of  Electrons  in  Silicon 

Dioxide,"  J.  C.  Ashley  and  V.  E.  Anderson,  IEEE  Trans.  Nucl. 

Sci.  NS-28,  4132-4136  (1981). 


B.  Summary  of  Mean  Free  Path  Calculations 


Results  for  electron  inelastic  mean  free  paths  in  ten  inorganic 
solids,  some  of  which  results  were  not  readily  available  in  the  open 
literature,  were  collected  and  presented  in  a  particularly  simple  form. 

Values  of  k  and  p  were  given  for  determining  the  mean  free  path  X  as  a 
function  of  electron  energy  E  from  the  simple  power-law  expression 
x  =  kEp  for  values  of  E  from  200  eV  through  10  keV.  The  resulting  table 
of  k-p  values  and  comparisons  of  mean  free  paths  with  experimental  data 
and  with  other  theoretical  results  were  given  in  the  following  paper: 

"Electron  Inelastic  Mean  Free  Paths  in  Several  Solids  for  200  eV  <  E  *  10  keV, 
J.  C.  Ashley  and  C.  J.  Tung,  Surf.  Interface  Anal.  4^,  52-55  (1982). 

C.  Polyethylene 

Previous  results,  derived  using  an  insulator  model,  for  the  energy 
loss  per  unit  path  length  and  inelastic  mean  free  path  of  an  electron 
with  an  energy  *10  keV  [ Radiat .  Res.  83,  1-18  (1980)]  were  modified  by 
including  an  exchange  correction,  a  sum-rule  constraint,  and  a  small 
relativistic  correction.  These  new  results  were  compared  with  the 
previous  results,  with  tabulated  values  of  stopping  power  based  on  Bethe 
theory,  and  with  predictions  from  our  universal  formula  for  organic 
solids  in  the  publication  "Energy  Losses  and  Inelastic  Mean  Free  Paths 
of  Low-Energy  Electrons  in  Polyethylene,"  J.  C.  Ashley,  Radiat.  Res.  90, 
433-436  (1982). 

D.  Condensed  Organic  Materials 

A  simple,  "nne-mode"  approximation  for  the  wavevector-  and  frequency- 
dependent  energy-loss  function  was  employed  to  derive  an  expression  for 


-5- 


the  electron  inverse  mean  free  path,  a*1,  including  the  effect  of  exchange 
between  the  incident  electron  and  the  electrons  in  the  medium.  The 
inverse  mean  free  path  in  a  specific  material  can  then  be  obtained  by 
averaging  a-1  over  the  optical  oscillator  strength  distribution  which  is 
proportional  to  u>Im[-l/£(u>)].  For  materials  for  which  the  optical 
oscillator  strength  distribution  may  be  approximated  by  a  single,  Drude- 
type  function,  especially  organic  materials,  the  general  results  were 
cast  in  a  particularly  simple  form.  With  four  constants  obtained  from 
optical  data  on  a  given  material,  electron  mean  free  paths  in  that 
material  can  be  predicted  for  electron  energies  *150  eV.  An  approximate 
expression  was  also  deduced  for  estimating  electron  mean  free  paths  in 
materials  for  which  no  optical  data  is  available.  Mean  free  paths 
determined  from  this  approximate  expression  agree  well  with  those 
obtained  from  our  earlier  universal  curve  predictions.  For  polyethylene, 
excellent  agreement  between  results  from  this  model  and  those  of  a  more 
detailed  insulator  model,  as  described  in  section  C,  gave  additional 
support  for  this  simple  theoretical  approach.  Agreement  between  theoretical 
predictions  and  experimental  measurements  of  electron  mean  free  paths  in 
organics,  which  in  some  cases  is  quite  good,  is  on  the  whole  less  than 
satisfactory.  These  results  are  described  in  detail  in:  "Simple  Model 
for  Electron  Inelastic  Mean  Free  Paths:  Application  to  Condensed  Organic 
Materials,"  J.  C.  Ashley,  >r.  Electron  Spectroec.  Relat.  Phenom.  28, 

177-194  (1982). 

E.  Additional  Model  Calculations 

We  have  continued  to  examine  simple  models  for  the  response  of 
solids  to  energy  and  momentum  transfers  from  low-energy  electrons 
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(slO  keV)  for  prediction  of  inelastic  mean  free  paths  (IMFP)  and  energy 
losses  per  unit  pathlength.  Some  of  these  results  for  the  IMFP  were 
discussed  in  section  1-0  of  this  report.  Here  we  discuss  predictions  of 
stopping  powers.  We  take 


Im[(-l/c(q,  w)]  = 


dui'w'  Im[-l/e(o, 


6^0)  -  (w"  +  q2/2)^/u> 


(1) 


o 

for  the  energy-lots  function  which  assumes  a  simple  quadratic  extension 
from  the  optical  limit  (q=0)  into  the  energy  and  momentum  transfer 
plane.  The  energy  loss  per  unit  pathlength,  in  atomic  units,  is  calc  -lued 
from 


Im[-l/e(q, 


to) 


(2) 


where  u>  and  q  for  w  >  0  are  restricted  to  values  for  which  w  <  q(v-q/2), 
and  v  is  the  electron  speed.  For  a  fixed  u>,  the  values  of  q  lie  in  the 
range  q_  s  q  s  q+  where 

q+  =  v  +  v^v2  -  2w  .  (3) 


Equations  (l)-(3)  lead  to 


Jr  f 

itV2  J 


du'u)'  Im[-l/e(o,  w')] 


F(f-, 


u>) 


where 


u>)  =  o(w 


qf 

■  T 


u)')6(u)' 


(4) 


(5) 


and  e(x)  *  1  for  x  >  0  and  g(x)  =  0  for  x  <  0.  We  define  an  exchange- 

corrected  differential  inverse  mean  free  path  t  by 

exc  * 
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/“ 

dw'w'  Im[-l/e(o,  w')]  |f(|— ,  w',  uj) 


+  F (5— > 


0 

2 

,  w 


(6) 


v2  ,  -  x  rF/v2  -  \F/v2  -  v2  ,  .  \-ii/2j. 

+  0)  -  U)J  -  L'  V'5 — 9  w  »  '9  w  9  Ty~~  +  w  •  w;J 


T 


2~» 


r 


where  the  expression  in  {  }  reduces  to  a  form  proportional  to  the  non- 

v2 

relativistic  Miller  cross  section  when  w  and  ^ —  u>  are  large  compared 
to  the  "binding  energy"  u'.  With  Eq.  (6),  S  is  given  by 


where  <*>_  s  J  J  D  +  “  f]  ~  3  • 

The  integrations  over  w  in  Eq.  (7)  can  be  done  analytically  resulting 
in  the  following  expression: 


1  y2 

s  =  ~2  f  dw'u>'  Im[-l/e(o,  «')]  Tn 

0 

(sin  1  •tn[(l+a/2  -  Va(l+a/4)  )/a]|  . 


(8) 


In  this  expression  a  =  2 <*>7v2,  s  =  (l-2a) 1/2,  and  ^  is  an  elliptic 
integral  of  the  first  kind  defined  by 


-8- 


k)  = 


sin* 

dx  [<1  -  x2)(l  -  k2x2)]'1/2 


(9) 


Given  a  set  of  optical  data  over  the  required  energy  range,  the 
integral  over  <d'  in  Eq.  (8)  could  be  done  numerically  to  obtain  S. 

However,  a  simple  expression  can  be  obtained  by  expanding  the  analytic 
functions  resulting  from  the  <d  integration  in  powers  of  cd'/v2.  If  we 
assume  v2/ 4  is  large  compared  to  the  value  of  id'  for  which  <D'Im[-l/e(o,w')] 
is  largest,  then  only  the  first  few  terms  in  the  expansion  should  be 
sufficient  to  approximate  S. 

Through  second  order  terms  in  id'/v2,  we  find 


1  v2 
'2  2 


*  f2  2  du'w'  Im[-l/e(o,  w')] 1[£m  |  ^ 


(10) 


ru>'  ,3/2  1  /  id  ' , 2  n  2v2  31  /id',2- 

( ln  ~  '  T7  ■ 


From  this  form,  in  the  large  v  limit,  we  obviously  recover  the  non- 
relat'vistic  Bethe  formula 


SB,NR 


(11) 


where  a*  =  4TrnQZ,  nQ  is  the  number  of  atoms  or  molecules  per  unit 
volume  with  Z  electrons  per  atom  or  molecule  and  I  is  the  mean  excitation 
energy  defined  by 


/ 


dio'io'  Im[-l/e(o,  u>')]  ■Znu 


//■- 


Im[-l/e(o,  cd')3  ,  (12) 


The  stopping  power  has  been  evaluated  from  Eq.  (10)  for  polyethylene 


using  our  optical  data  [ Radiat .  Res.  83,  1-18  (1980)],  including  carbon 
K-shell  contributions,  with  the  results  shown  in  Table  I  ("this  result"). 


Table  I.  Stopping  Power  of  Polyethylene 


E(eV)  _ 

0 

Stopping  Power  in  eV/A 

This  Result 

"R.R." 

"Bethe" 

80 

2.68 

3.00 

2.12 

100 

2.92 

3.14 

2.67 

150 

2.92 

3.05 

2.85 

200 

2.72 

2.82 

2.73 

300 

2.30 

2.37 

2.38 

400 

1.97 

2.05 

2.08 

600 

1.56 

1.64 

1.66 

800 

1.31 

1.38 

1.40 

1000 

1.14 

1.20 

1.21 

1500 

0.882 

0.917 

0.917 

2000 

0.724 

0.757 

0.747 

4000 

0.437 

0.454 

0.445 

6000 

0.320 

0.331 

0.324 

8000 

0.255 

0.263 

0.258 

10000 

0.213 

0.221 

0.216 

These  results  are  in  good  agreement  with  our  recent  results.  Column 
"R.R.,"  Table  I,  obtained  using  an  insulator  model  to  describe  the 
valence  electrons  and  generalized  oscillator  strengths  for  K-shell 
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ionization  (see  section  C).  Surprisingly,  both  sets  of  results  agree 
with  the  Bethe-theory  predictions,  Eq.  (11),  (Column  "Bethe,"  Table  I), 
to  within  a  few  percent  down  to  ^200  eV.  Part  of  the  difference  in  the 
S  values  in  Table  I  between  "this  result"  and  "R.R."  is  due  to  a  dif¬ 
ference  in  I  in  the  two  cases.  From  the  energy-loss  function  obtained 
from  optical  data,  I  =  62  eV.  The  insulator-model  results  ("R.R.") 
involved  a  fit  of  a  sum  of  terms  to  Im[c(o,6j)].  This  function,  plus  the 
Kramers-Kronig  relation  between  the  real  and  imaginary  parts  of  e,  were 
used  to  determine  the  valence-electron  contribution  to  Im[-l/e(o,u>)]. 

This,  plus  the  K-shell  contribution  as  determined  from  generalized 
oscillator  strengths,  resulted  in  I  =  58  eV.  This  slightly  smaller  I 
corresponds  to  slightly  larger  values  of  S  than  for  I  =  62  eV,  e.g.,  a  7% 
difference  in  I  gives  about  a  1.5%  difference  in  S  at  8000  eV  according 
to  Eq.  (11). 

The  advantage  of  this  model  calculation  is  that  a  reasonably  simple 
function,  when  integrated  over  the  optical  energy  loss  function,  yields 
S  from  Eq.  (10).  The  disadvantage  is  that  the  optical  energy-loss  function 
is  not  available  for  all  materials  of  interest.  A  more  general  problem 
is  that  almost  no  measurements  are  available  of  stopping  powers  for 
electrons  with  energies  less  than  10  keV  as  a  check  on  the  theoretical 
descriptions. 


II.  ELECTRON  ATTENUATION  LENGTHS  IN  CARBON 


Measurements  have  been  made  of  electron  attenuation  lengths  as  a 
function  of  incident  electron  energy  over  the  energy  range  from  100  to 
3000  eV  for  electrons  in  amorphous  carbon  films.  These  measurements 
were  undertaken  as  a  direct  check  on  the  values  predicted  by  our  universal 
expression  for  the  mean  free  path  of  electrons  in  solid  organic  insulators 
[. Badiat .  Res.  81_,  364-373  (1980);  IEEE  Trane.  Hud.  Sai.  NS-27,  1454-1458 
(1980)]. 

The  attenuation  lengths  were  obtained  from  electron  transmission 
measurements  through  self-supporting  carbon  films.  The  films  were 
prepared  by  the  arc  evaporation  of  high  purity  carbon  electrodes  in 
vacuum.  Carbon  was  deposited  onto  glass  slides,  which  had  been  spin- 

o 

coated  with  a  dilute  soap  solution,  at  a  rate  of  approximately  5  A/sec. 

With  the  soap  acting  as  a  parting  agent,  the  carbon  films  were  floated 
off  the  glass  slides  onto  distilled  water,  and  then  picked  up  on  829- 
line-per-inch  Ni  mesh.  The  films,  which  had  thicknesses  in  the  range 

O 

from  60  to  200  A,  were  examined  using  an  optical  microscope  to  locate 
continuous  segments  which  might  be  suitable  for  the  electron  transmission 
measurements.  Since  electrons  could  be  transmitted  freely  through  holes 
too  small  to  be  resolved  optically,  only  those  segments  for  which  the 
electron  transmission  became  unmeasurable  at  low  incident  electron 
energies  were  used  in  the  measurement  of  electron  attenuation  lengths. 

Electrons,  accelerated  by  an  electron  gun,  were  directed  at  the 
sample,  and  the  electrons  transmitted  in  a  relatively  large  solid  angle 
about  the  forward  beam  direction  were  collected.  In  order  to  discriminate 
against  the  majority  of  the  electrons  which  had  been  scattered  inelastically. 
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the  detector  incorporated  a  retarding  potential  grid  which  was  set  at 
ten  volts  less  than  the  electron  gun  accelerating  potential.  For  each 
incident  electron  energy,  the  e^ectror  attenuation  it-nct.n,  L,  associated 
with  inelastic  scattering  was  then  obtained  from  the  electron  transmittance, 
T,  using  the  relationship 


T  -  e-t/L  _  1  I  *BG 
c  I  -  I 

lo  1BG 


where  t  is  the  sample  thickness,  I  is  the  incident  electron  intensity, 

I  is  the  transmitted  electron  intensity,  and  IBG  is  the  background 
intensity. 

The  experiments  were  performed  in  an  ion-pumped,  ultra-high  vacuum 
chamber.  To  reduce  the  effect  of  magnetic  fields  on  the  electron  beam, 
a  mu-metal  shield  enclosed  the  sample  region.  Holes  were  cut  in  the 
shield  for  the  electron  gun,  the  detector  and  target  assemblies,  and  for 
an  observation  window.  Great  care  had  to  be  used  in  the  initial  pump- 
down  to  avoid  rupturing  the  fragile  films.  An  oil -free  mechanical  pump, 
pumping  through  a  small  tube,  was  used  first,  followed  by  sorption  and  ion 
pumps.  The  ultimate  vacuum  was  about  10-8  Torr,  increasing  to  5xlo-8 
Torr  when  the  electron  gun  was  in  operation. 

The  thicknesses  of  the  films  were  determined  from  the  measured 
transmi ttance  of  He-Ne  laser  photons,  using  a  transmittance  versus 
thickness  curve  which  had  been  calibrated  using  a  quartz  crystal  thickness 
monitor.  The  thicknesses  were  76,  98,  103,  and  163  A  with  experimental 
uncertainties  in  these  numbers  of  ±20%. 

The  electron  attenuation  lengths  obtained  for  the  four  carbon  films 
used  in  this  study  are  shown  in  Fig.  1 1-1  as  a  function  of  electron 
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energy.  Representative  error  bars  are  shown,  and  these  result  primarily 
from  the  20%  uncertainty  in  the  film  thicknesses.  The  dotted  curves  are 
calculated  using  our  "universal"  expression  for  electron  mean  free  paths 
in  solid  organic  insulators  [IEEE  Trans.  Nunl.  Sci.  MS-27,  1454-1458  (1980)] 
for  carbon  film  densities  of  1.8,  2.0,  and  2.2  g*cm-3.  The  dashed  curves 
are  calculated  using  Penn's  equation  [J.  Electron.  Spectrosc.  Relat. 

Phenom.  9,  29-40  (1976)]  based  on  a  free  electron  gas  model. 

To  make  a  valid  comparison  between  experiment  and  theory  it  is 
necessary  to  know  the  density  of  the  amorphous  carbon  films  used  in  the 
electron  transmission  measurements.  A  preliminary  determination  of  this 
density  gave  2.1  g*cm-3.  Assuming  this  value  of  density  to  be  correct, 
examination  of  Fig.  1 1- 1  shows  that  the  experimental  values  of  electron 
attenuation  length  are  in  good  agreement  with  the  predictions  of  electron 
mean  free  path  obtained  from  our  universal  curve  from  ^500  to  3000  eV. 

Below  ^500  eV,  the  experimental  values  of  electron  attenuation  length 
are  larger  than  predicted  from  the  universal  curve.  Furthermore,  the 
observed  dependence  of  electron  attenuation  length  on  energy  is  not  that 
which  is  predicted.  One  reason  for  the  differences  between  experiment 
and  theory  for  incident  electron  energies  $500  eV  may  be  that  some 
inelastically  scattered  electrons  may  be  included  in  the  measured 
transmitted  intensity,  I.  Alternatively,  the  universal  curve  may  not 
apply  in  this  energy  region.  Further  studies  are  proposed  to  clarify 
the  questions  raised  by  this  study. 


III.  CHARGED-PARTICLE  TRACK  STRUCTURE  IN  NON-METALLIC  SOLIDS 


Our  core  plasma  model  of  the  expanding,  partially  compensated  posi¬ 
tive  ion  track  created  when  a  densely  ionizing  charged  particle  proceeds 
through  matter  has  been  generalized  so  that  the  collective  effects  are 
included  a  priori.  The  time  evolution  of  delta  rays  that  are  generated 
in  the  primary  physical  stage  of  the  energy  deposition  process  is 
followed  by  means  of  an  age- type  approximation  that  includes  the  self- 
Coulomb  forces  of  the  delta-ray  distribution  as  well  as  the  Coulomb 
attraction  of  the  positive  ion  column.  On  a  slower  time  scale  the  ion 
column  expansion  is  described  self-consistently  in  a  hydrodynamical 
framework  and  the  atomic  recoil  distribution  generated  in  the  expansion 
is  modeled  on  a  multi-group  diffusion  approximation.  Calculations  with 
this  scheme  will  furnish  a  complete  description  of  the  evolution  in  time 
and  space  of  the  energy  distribution  of  atomic  recoils  and  the  ultimate 
conversion  of  this  kinetic  energy  into  thermal  energy  in  a  detailed 
space-time  representation.  Monte  Carlo  checks  of  these  calculations  are 
planned.  Some  results  of  this  work  have  been  published  in:  "A  Core 
Plasma  Model  of  Charged  Particle  Track  Formation  in  Insulators," 

R.  H.  Ritchie  and  C.  Claussen,  Nuol.  Instrun.  Methods  198,  133-138 
(1982). 
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IV.  ELECTRON  INTERACTIONS  IN  Si  and  Si02  IRRADIATED  BY  PHOTONS 

Preliminary  calculations  were  carried  out  of  the  relative  energy 
densities  deposited  as  a  luiKtion  oi  depth  by  monoenergetic  x-ray  photons 
incident  on  a  0.5-pm  slab  of  Si02  sandwiched  between  two  slabs  of  Si. 

The  slabs  were  infinite  in  lateral  extent  and  the  photons  entered  as  a 
parallel  beam  perpendicular  to  the  planes  of  the  Si-Si02-Si  interfaces. 

The  front  Si  layer  was  3-pm  thick  and  the  back  layer  was  essentially 
infinitely  thick.  Photon  energies  of  8.04  keV,  17.4  keV,  and  60.0  keV 
were  used.  The  first  two  energies  are,  respectively,  those  of  the 
characteristic  x  rays  of  Cu  and  Mo.  The  photoelectric  effect  dominates 
in  these  two  cases.  The  energy  60.0  keV  was  chosen  as  an  example  in 
which  photoelectric  absorption  and  Compton  scattering  are  about  equally 
probable  both  in  the  case  of  Si  and  of  Si 02 . 

Figure  I V- 1  shows  the  energy  deposited  per  unit  volume  (arbitrary 
units)  in  the  array  as  a  function  of  depth  for  irradiation  by  8.04  keV 
photons  incident  from  the  left.  The  Si02  slice  is  located  between  the 
depths  of  3.0  pm  and  3.5  pm.  Rather  than  selecting  the  photon  interaction 
sites  at  random,  we  let  them  occur  at  uniform  intervals  in  the  Si  and 
Si02,  the  linear  density  being  adjusted  in  the  ratio  of  the  total  macro¬ 
scopic  cross  sections  for  photon  interaction.  Since  our  primary  purpose 
is  to  study  the  effects  of  the  interface  between  the  two  materials,  it 
Is  not  necessary  to  have  photon  interactions  at  all  depths  in  the  array. 

To  calculate  the  data  for  Fig.  IV-1,  we  allowed  photons  to  interact  only 
between  the  depths  of  about  2.5  pm  and  4.0  pm,  indicated  by  the  vertical 
arrows  on  the  abscissa  in  Fig.  IV-1.  These  depths  were  sufficiently  far 
from  the  Si02  layer  to  establish  electronic  equilibrium  in  the  Si  on 
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bo  th  sides  of  it.  The  results  show  that  the  energy  deposited  per  unit 
volume  drops  as  one  goes  from  the  Si  into  the  Si02  layer  in  the  direction 
the  photons  are  traveling.  The  energy  deposited  per  unit  volume  then 
rises  in  going  from  the  Si02  into  the  Si  on  the  back  side.  These  results 
show  a  gradient  in  energy  deposition  across  the  Si -Si 02  interface.  A 
total  of  4000  photon  interactions  were  used  in  obtaining  the  data  for 
Fig.  IV-1. 

Figure  IV-2  shows  preliminary  results  for  the  17.4  keV  x  rays.  In 
this  case  1800  photon  interactions  were  used.  An  energy  deposition 
gradient  in  the  Si02  layer  is  also  seen  here. 

Figure  IV-3  shows  the  raw  data  calculated  for  the  60-keV  photons. 

The  statistics  are  poor,  as  only  450  interactions  were  used  spread  over 
the  interval  from  0.1  to  6.4-pm  depth.  If  one  combines  the  data  for  the 
Si02  in  Fig.  IV-3  in  intervals  of  0.2  pm  instead  of  the  0.1  pm  shown, 
there  is  evidence  of  a  gradient  similar  to  those  found  at  the  two  lower 
photon  energies. 

Based  on  these  figures,  final  decisions  will  be  made  on  the  geometrical 
configurations  and  photon  energies  of  most  interest.  Final  computer 
runs  will  then  be  made  with  a  sufficient  number  of  photon  interactions 
to  assure  satisfactory  statistical  resolution  in  the  results. 
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